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ABSTRACT 

We present far-UV observations of the Crab nebula and its pulsar made with the Space Telescope 
Imaging Spectrograph onboard the Hubble Space Telescope. Broad, blueshifted absorption arising in the 
nebula is seen in C IV A1550, reaching a blueward velocity of ~ 2500 km s . This can be interpreted as 
evidence for a fast outer shell surrounding the Crab nebula, and we adopt a spherically symmetric model 
to constrain the properties of such a shell. From the line profile we find that the density appears to 
decrease outward in the shell. A likely lower limit to the shell mass is ~ 0.3 Mq with an accompanying 
kinetic energy of ~ 1.5 x 10 49 ergs. A fast massive shell with 10 51 ergs cannot be excluded, but is less 
likely if the density profile is much steeper than p(R) oc R~ 4 and the maximum velocity is <J 6000 km s — 1 . 

The observations cover the region 1140— 1720 A, which is further into the ultraviolet than has previously 
been obtained for the pulsar. With the time-tag mode of the spectrograph we obtain the pulse profile in 
this spectral regime. The profile is similar to that previously obtained by us in the near-UV, although the 
primary peak is marginally narrower. Together with the near-UV data, and new optical data from the 
Nordic Optical Telescope, our spectrum of the Crab pulsar covers the entire region from 1140- 9250 A. 
Dereddening the spectrum with a standard extinction curve we achieve a flat spectrum for the reddening 
parameters E(B — V) = 0.52, R = 3.1. This dereddened spectrum of the Crab pulsar can be fitted by a 
power law with spectral index ct v = 0.11 ± 0.04. The main uncertainty in determining the spectral index 
is the amount and characteristics of the interstellar reddening, and we have investigated the dependence 
of a u on E(B — V) and R. In the extended emission covered by our 25" x 0'.'5 slit in the far-UV, we 
detect C IV A1550 and He II A1640 emission lines from the Crab nebula. Several interstellar absorption 
lines are detected along the line of sight to the pulsar. The Lya absorption indicates a column density 
of (3.0 ± 0.5) x 10 21 cm~ 2 of neutral hydrogen, which agrees well with our estimate of E(B — U)=0.52 
mag. Other lines show no evidence of severe depletion of metals in atomic gas. 

Subject headings: pulsars: individual (Crab pulsar) — ultraviolet: stars — ultraviolet: ISM — dust: 
extinction — supernova remnants — instrumentation: spectrographs 



1. INTRODUCTION 



The Crab nebula and its pulsar (PSR 0531+21) are 



among the mo3t studied objects in the sky. The di acovc ry - 
jf the Crab pulnai' an a taut rotating radi o pulociiP ( |Staiolin | 



Roifonetoin 1968 -- Comolla ot al. 1969| ) paved the 
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light: the Cr ab pulsar (|Cocke et al. 196S|), the LMC pul- 
sar 054 0-69 QMiddleditch fc Pcnnypackcr 1985|), the Vela 
pulsa r ( IWallace et al. 19771 ), p SR 0656+14 (IShcarcr et"aT 



1997) and (possibly) the Geminga pulsar (Shearer et al 



1998). In the near-UV (NUV), pulsations have only been 



Also, the position of the Crab pulsar in the con 



Gull 



et al. 1998 , henceforth G98). Due to the faintness of these 



ter of the Crab nebula, which is the remnant of supernova 
1054, clearly supports the supernova - neutron star con- 
nection. Soon after the radio detection the pulsar was also 
shown to emit optical pulsations (Cocke, Disney, & Taylor 
1969). This established that the pulsating star was the 
well known south preceding star in the center of the neb- 
ula, which early opti cal spectroscopy showed to emit a fea- 
tureless continuum (Minkowski 1942). To date, more than 
1000 radio pulsars are known, but only the following few 
have optical counterparts known to pulsate also in visible 



objects in the optical and in the ultraviolet, the spectro- 
scopic information is very limited. PSR 0540-69 was ob- 
served with the Faint Object Spectrograph (FOS) onboard 
Hubble Space T elescope (HST) in the 2500 - 5000 A range 
( Hill et al. 1997 ) and showed a rather steep power law spec- 
trum. These observations were, however, contaminated 
by nebular emission. The Geminga pulsar was observed 
with the Keck telescope (Martin, Halpern, & Schiminovich 
1998), but the spectrum has very low signal-to- noise be- 
cause this pulsar is exceedingly faint in the optical. The 
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only pulsar for which good signal-to-noise spectroscopy in 
the optical and ultraviolet can be obtained is the Crab 
pulsar. Surprisingly enough, very little has been done in 
this respect since the optical observations of Oke (1969). 
In particular, until the study by G98, no UV spectroscopy 
of the Crab pulsar had been published since the first at- 
tempts by the Intern ational Ultraviolet Explorer (IUE) 
dBcnvenuti ct al. 1980| ). The IUE data cover only the NUV 
region (2000 — 3150 A) and have poor signal-to-noise. The 
HST/STIS (Space Telescope Imaging Spectrograph) data 
from G98, and the new data presented here, clearly super- 
sede these early attempts. 

The Crab pulsar has been extensively studied over a 
very broad wavelength range, from the radio up to 7- 
rays (e.g., Lyne fc Graham-Smith 1998 ). The high en- 



ergy emission, from infrared (IR) to 7-rays, is believed 
to be the result of the same emission mechanism (e.g., 
Lyne fc Graham-Smith 1998 ). It is therefore of interest 
to fill in the gaps in the observed spectrum of the pulsar 
in this range. Although the pulsar is relatively bright in 
the optical, UV observations are difficult due to the large 
extinction toward the C rab, E(B — V) ~ 0.5 mag (e.g., 
Davidson fc Fesen 1985| , and references therein). Here, 
we present UV observations of the Crab pulsar further 
into the far-UV (FUV) (1140 - 1720 A) than has previ- 
ously been obtained. These are presented together with 
our previous NUV-data (1600 - 3200 A)(G98) and new 
optical data from the Nordic Optical Telescope (NOT). 
Due to the large extinction correction, great care must be 
taken to draw conclusions about the intrinsic spectrum, 
and thus the emission mechanism of the pulsar. However, 
this procedure might also give a hint on the absorption 
properties of the dust in the direction toward the pulsar. 

In addition to the pulsar emission, we detect emission 
lines from the Crab nebula itself in the FUV. In particular, 
the strength of the C IV A1550 emission can be of inter- 
est for abundance determinations. Even more interesting 
is the broad C IV A1550 absorption line from the nebula 
detected against the pulsar continuum. This line provides 
information on the nature of the SN 1054 event. 

Although the Crab nebula has been studied extensively, 
the nature of the progenitor remains unknown. According 
to models based on the existence of the central neutron 
star, as well as on nucleosynthesis arguments, the zero- 
age main sequence (ZAMS) mass of the progenitor was 
probably in the range 8 — 13 Mq (Nomoto 1985). The 
amount of material observed in the nebula (4.6 ±1.8 Mq) 
seems too low to account for this (Fesen, Shull, & Hur- 
ford 1997). Furthermore, the velocities of the filaments 
(~ 1400 km s -1 , Davidson & Fesen 1985) give an uncom- 
fortably low kinetic energy (<; 1 x 10 50 ergs) compared to 
other supernova remnants, i.e., at least an order of mag- 
nitude less than the canonical energy of supernovae, 10 51 
ergs. 

The "missing mass" could either be in a slow progen- 
itor wind, or in a fast, hitherto undetected, shell ejected 
at the explosion (Chevalier 1977). If the latter is true, as 
is hinted by the observations of the outer [O III] skin of 
the nebula (Hester et al. 1996; Sankrit & Hester 1997), 
this shell might account for the missing mass and kinetic 
energy of the nebula. The question remains, however, why 
such a shell has escaped detection despite many efforts to 



observe it (see, e.g., Fesen et al. 1997 ). One possibility 
is that the low density of the surrounding gas is not high 
enough to give rise to detectable circumstellar emission 
when interacting with the ejecta; neither X-ray nor ra- 
dio searches have indicated any evidence of circumstellar 
interaction between fast ejecta and ambient gas (Mauche 
& Gorenstein 1989; Frail et al. 1995). If a fast shell is 
absent, the birth of the Crab was definitely a low energy 
event. This would call for a revision of our understanding 
of supernova explosions, especially since SN 1054 appar- 
ently was not unusually dim according to historical records 
(Chevalier 1977; Wheeler 1978). 

It is thus of great interest to further investigate 
whether there is a stellar wind or supernova ejecta out- 
side the observed nebula, and what velocity this gas may 
have. Lundqvist, Fransson, & Chevalier (1986, henceforth 
LFC86) proposed to search for a fast shell by looking in the 
UV toward the Crab pulsar. Their time dependent pho- 
toionization calculations showed that C IV A1550 could 
show up in blueshifted absorption if the ionization history 
of the shell was as predicted in some models of Reynolds 
& Chevalier (1984). Here, we present the detection of this 
broad absorption line, and discuss its implications for the 
fast shell around the Crab nebula. 

First we discuss the observations and reductions (§2). 
We then (§3) discuss the pulse profile, the amount of red- 
dening toward the pulsar and the intrinsic pulsar spectrum 
in the optical/UV. In §3 we also discuss the lines originat- 
ing from the interstellar gas toward the pulsar and from 
the Crab nebula itself. Some of these observations con- 
strain the properties of a possible outer shell. In §4 we 
summarize our conclusions. 

2. OBSERVATIONS, REDUCTIONS AND RESULTS 

2.1. HST Far-UV Observations 

The Crab pul sar was observed on January 22, 1999 us- 
ing HST/STIS ( |Kimble et al. 1998|) with the FUV Multi 
Anode Micro-channel Array (MAMA) detector. The low 
resolution grating G140L was used, which covers the wave- 
length interval 1140 — 1720 A. These observations were 
made in the time-tag mode and used a slit of 52" x Of! 5. 
The spectral resolution is 0.58 A pixel -1 , and the plate 
scale is 0'.'0244 pixel" 1 . This means that only 25" of the 
long-slit is actually projected onto the detector. In to- 
tal, six orbits of observations, including target acquisition, 
were used. These were divided into two visits. A log of the 
observations is shown in Table 1. The total on-target ex- 
posure time amounted to 14,040 seconds. The orientation 
of the slit is shown in Figure 1. 

2.1.1. Time-resolved emission 

The time-tag mode on the STIS allows us to resolve 
the emission from the Crab pulsar both in wavelength and 
time. The time resolution obtained in this mode is 125 /is. 
As these observations were the first to utilize the time- 
tag capabilities of HST/STIS in the FUV for a known 
periodic variable, special software, developed at Goddard 
Space Flight Center (GSFC) was used to obtain the pulse 
profile for the pulsar. The analysis followed the procedures 
outlined in G98. For each of the six datasets, a time aver- 
aged image was produced to trace the position of the pul- 
sar spectrum. A 13 pixel wide window was used to extract 
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events in the pulsar spectrum as well as in the background 
emission at both sides of the pulsar emission. The arrival 
time of each 125 //s sample was converted to a solar sys- 
tem barycenter arrival time. The position of the HST with 
respect to earth center was computed by the Flight Dy- 
namics Facility at GSFC with errors less than 200 meters. 
The position of the earth with respect to the solar system 
barycenter was computed using a routine, SOLSYS, sup- 
plied by the U. S. Naval Observatory (Kaplan et al. 1989). 
All events in the pulsar spectrum (and in the background 
regions) were assigned a (barycentric) arrival time with re- 
spect to the start of the first exposure. As we found a small 
drift in the times recorded in the FITS headers, we used 
the internal clock from the engineering mode header to do 
this. To determine the period we folded the arrival times 
modulo a grid of test frequencies, and corrected for the 
slowdown rate of the pulsar. Pulse profiles were calculated 
as histograms of the function f(t) = fit + ft 2 /2, where the 
data were coadded into 512 phase bins. The appropriate 
value of /, was then determined by maximizing the sum 
of squares of the values in the pulse profile. 

In this procedure we used a value for / from ra- 
dio observations at Jodrell bank (Lyne, Pritchard, & 
Roberts 1999), while the second time derivative is unim- 
portant for this purpose. This resulted in a measured pe- 
riod of P=33.492675 ms at Modified Julian Date (MJD) 
51200.549, the time of the beginning of our first obser- 
vation (number O4ZP01010). As the data were obtained 
during a time period of eight hours, we could also deter- 
mine the pulsar slowdown rate from our observations. We 
obtained P = (4.0 ±0.4) x 10~ 13 s s" 1 , which is consistent 
with the value used above from the radio observations, 
4.2 x 10~ 13 s s" 1 . 

In Figure 2 we show the Crab pulsar pulse profile in the 
FUV regime. It was obtained by subtracting the back- 
ground from the pulse profile obtained for the period given 
above. For comparison, the figure also shows the NUV 
pulse profile from G98. 

2.1.2. Phase- averaged spectrum 

Averaging over the pulse we obtain a phase-averaged 
spectrum of the pulsar which covers the region 1140 — 
1720 A. This spectrum probes the emission of the Crab 
pulsar further into the UV than has previously been done 
and is shown in Figure 3, together with the NUV spec- 
trum of G98. These spectra overlap nicely and cover to- 
gether the whole range from 1140 — 3200 A. The combined 
spectrum offers the possibility to deduce the amount and 
characteristics of the interstellar reddening, as well as to 
determine the dereddened pulsar spectrum itself. 

The FUV-spcctrum was extracted with a 13 pixel wide 
window. The reductions were made using the CALSTIS 
software developed at the GSFC. These IDL-routines flat- 
field the images and then the point source spectrum is 
localized and traced on the detector. The extracted pul- 
sar spectrum is background subtracted and converted to 
absolute flux units using the G140L sensitivity table. The 
accuracy of the absolute flux calibration is ~ 15% over 
the full wavelength scale. Wavelengths are assigned from 
a library dispersion solution, while the zero point adjust- 
ments are determined from arc frames taken through the 
0'.'05 slit for each science observation. The wavelengths are 



then converted to heliocentric wavelengths. The accuracy 
of the wavelength solution is about 0.4 A. 

2.2. Optical observations 

In addition to the UV spectrum, we have also collected 
data in the optical regime. During several nights in De- 
cember 1998 we did spectroscopy of the Crab pulsar us- 
ing the Andalucia Focal Reducer and Spectrograph (AL- 
FOSC) at the 2.56m NOT on La Palma. In total we ob- 
tained 11.25 hours of data in five different grisms. The l'/2 
slit was used for all observations (see Table 2 for more de- 
tails) . Not all nights were photometric and the seeing was 
generally just above 1". The data were bias subtracted 
and flatfielded. Wavelength calibrations were done using 
arc frames obtained with a helium lamp. Flux calibra- 
tion of the spectra was accomplished by comparison to 
the spectrophotometry: standard stars Feige 34 and G191- 
B2B. To avoid systematic errors due to background sub- 
traction the slit was put at two different position angles. 
All observations were made at low air masses and close to 
the parallactic angle to reduce the effects of atmospheric 
dispersion; the Crab pulsar passes just 7° from zenith as 
viewed from La Palma in December. The slit positions for 
the NOT observations are shown in Figure 1. 

2.3. Combined UV / optical spectrum 

The combined optical and UV observations (both NUV 
and FUV) are shown in Figure 4. It covers the region 
1140 — 9250 A. Although great care was given to the back- 
ground subtraction of the nebula, the optical spectrum of 
the Crab pulsar was contaminated by over- and under- 
subtractions of strong nebular emission lines. These had 
to be taken out by hand, and we used the IRAF task 
SPLOT to interactively clean the spectrum. Points that 
deviated more than 4cr from a smooth continuum fit were 
rejected. This procedure was robust enough to exclude 
only clear cases of nebular contamination. The optical 
spectrum used is a combination of all the different spec- 
tra in Table 2. Note that an absolute flux calibration was 
not applied to the optical spectrum, as the observing con- 
ditions were often non-photometric. Instead we have ap- 
plied a grey shift to the spectrum to match the V-band 
observations of Percival et al. (1993) and Nasuti et al. 
(1996). As seen in Figure 4, this is well matched with the 
UV spectrum from the HST. 

3. DISCUSSION 
3.1. Pulse profile 

The work of Percival et al. (1993) showed small differ- 
ences in the optical versus NUV pulse profile shapes. They 
observed the Crab pulsar with the High Speed Photome- 
ter (HSP) onboard HST and found that the main pulse is 
slightly narrower in the UV than in the optical. Eikcn- 
bcrry et al. (1997) extended the analysis into the near-IR 
and found that the trend for a decreasing Full Width Half 
Maximum (FWHM) with decreasing wavelength seems to 
hold over the full UV-IR range. 

The time-tag mode of STIS/FUV-MAMA allows us to 
examine if the pulse profile of the emission is different in 
this wavelength region from that in the NUV. The most 
striking feature of the FUV pulse profile shown in Figure 2 
is certainly that it is very similar to the profile previously 
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obtained in the NUV (G98). It appears, however, that 
the primary peak is slightly narrower in the FUV than 
in the NUV, as indicated by the blow-up of that region 
in Figure 2. We measured the FWHM of the FUV and 
NUV primary peak to be 0.0405 and 0.0426 periods, re- 
spectively. The position of the peak was determined by 
a polynomial fit to the central 10 phase bins. The phase 
at half-maximum was then simply determined by linear 
interpolation between the two closest phase bins. To esti- 
mate a statistical error on the procedure used to measure 
the FWHM, we computed the FWHM for each of our six 
FUV observations. The standard deviation obtained in 
this procedure was 0.001 phase bins. The measured differ- 
ence in the FWHM of the primary peak can therefore be 
considered marginally significant. The secondary peak in 
Figure 2 might actually appear broader in the FUV. It is, 
however, much noisier than the primary peak and the same 
procedure as above could not determine any significant dif- 
ference to the 13% (3<r) level. The above findings are in 
agreement with the trend seen in Percival et al. (1993) 
and Eikenberry et al. (1997). The pulse period obtained 
from the FUV observat ions is P=33. 492 675 ms. Radio 
data from Jodrell Bank ( Lyne et al. 1999) ) determined the 
pulse period to be P=33. 492402 ms on January 15 1999. 
Using their values for the pulse period and its first time 
derivative on this date we can calculate the period at the 
time of our HST observations. The result is P=33.492676 
ms. This agrees with our estimate to 7 significant digits. 
The limiting errors in our computation of the period are 
the accuracy of the SOLSYS routine which has a barycen- 
tric velocity error of less than 8.0 x 10" 7 AU/day and the 
unknown accuracy of the rate of the STIS onboard clock. 

3.2. The UV extinction curve 

The phase averaged UV spectrum of the Crab pulsar 
must be corrected for a substantial amount of interstel- 
lar reddening. The value for E(B — V) has been es- 
timated by a number of authors; Wu (1981) obtained 
E(B -V) = 0.50 ± 0.03 mag by using the 2200 A dust ab- 
sorption feature, the nebular synchrotron continuum and 
an extinction curve derived from eight stars. Blair et al. 
(1992) also used the best fit to the UV nebular contin- 
uum and obtained E(B - V) = 0.51±g;g| mag. A differ- 
ent method was used by Miller (1973), who determined 
the reddening of the Crab nebula from observations of [S 
II] lines. Using modern values for the atomic parameters 
(Keenan et al. 1993; Ramsbottom, Bell, & Stafford 1996), 
and the extinction curve of Fitzpatrick (1999), his mea- 
surements gives E(B — V) = 0.50+0 06 ma 8- 

Our data of the pulsar itself allows us to estimate the 
value for E{B - V) by "ironing out" the 2200 A bump. To 
do this we assumed a standard value R — 3.1 and dered- 
dened the UV spectrum for different values of E(B — V). 
The dereddened spectra were then fitted by a power law 
and we chose the value for E(B — V) that minimized a av , 
the standard deviation of the power law fit in the region 
log ^=[14.98, 15.41], where the region including the Lya 
absorption was excluded. This procedure used the galac- 
tic mean extinction curve from Fitzpatrick (1999) and gave 
E(B — V) = 0.52 mag. This is in excellent agreement with 
the previous results stated above. As our data have better 
signal-to-noise and sampling than previous continuum fits, 



we will use this value as the best estimate of the extinction 
toward the Crab nebula throughout this paper. 

From the NUV spectrum taken of the Crab pulsar by 
IUE, claims were mad e for a peculiar extinction curve 
(Benvenuti et al. 1980). In particular, the 2200 A bump 
was reported to be substantially narrower than for the 
galactic mean extinction curve, a finding that could indi- 
cate that the supernova event itself had altered the grain 
composition in the Crab nebula. In principle both the 
extinction curve and the intrinsic pulsar spectrum are un- 
known, which of course makes it troublesome to disentan- 
gle these quantities. We take the following approach to 
this proble m: theoretical models favor a power law spec- 
trum (e.g., Lyne fc Graham-Smith 1998| ) and dereddening 
with E(B - V) = 0.52, R = 3.1 indeed gives a power law 
pulsar spectrum (Fig. 4), so we will simply assume the in- 
trinsic spectrum of the Crab pulsar to follow a power law 
F v = K(v/vo) a " ergs s" 1 cm" 2 Hz -1 . Here v is the fre- 
quency of the radiation and K is a constant that is nearly 
independent of a v when Uq is the logarithmic m ean fre- 
quency of the fitted bandpass ( Percival et al. 1993| ) . Using 
the extinction parameters above, E(B — V) — 0.52 and 
R = 3.1, the spectral index is a v = —0.035 in the UV 
range. 

By assuming that the intrinsic pulsar spectrum is indeed 
well represented by the obtained power law we derived an 
extinction curve toward the pulsar. This is shown in Fig- 
ure 7 together with the mean galactic extinction curve for 
R = 3.1 from Fitzpatrick (1999). The derived extinction 
curve has a moderately narrower dip (15%) and a some- 
what shallower rise in the extreme FUV. Apart from this, 
it overlaps nicely with the galactic mean extinction curve. 
Considering the large variety of measured UV extinction 
curves (see Fig. 2 in Fitzpatrick 1999), the derived ex- 
tinction curve toward the Crab can hardly be claimed to 
be peculiar. Assuming that the pulsar spectrum follows 
a power law, we conclude that we find no evidence for a 
non-standard extinction curve toward the Crab pulsar. 

3.3. The spectral index of the pulsar continuum 

According to models, the high energy pulsar emission, 
from IR to 7-ra ys, is produced by (curvature) s ynchrotron 

, and refer- 
ectrons fol- 
where E is 



Graham-Smith 1998 



radiation (e.g., Lyne 

ences therein). A number distribution of c 
lowing a power law N(E)dE = CE~~<dE, 
the energy, C is a constant and 7 is the electron spec- 
tral index, will produce synchrotron radiation that has a 
power-law distribution in flux density F v — K(y /vq) 01 " 
ergs s" 1 cm" 2 Hz" 1 . The photon spectral index, a„, is re- 
lated to the electron spectral index, 7, via a u = — (7— 1)/2 
for synchrotron radiation. 

The early low resolution spectroscopy of Oke (1969) 
appears to peak in the middle of the observed region 
(~ 3400 - 8000 A). He reports a slope of a v = -0.2 with 
no stated errors, although he cautions that the uncertainty 
in the reddening correction could allow also for a positive 
slope. Much of the theoretical work on the optical emission 
mechanism of pulsars has been b ased on this finding (cf. 
Ginzburg & Zheleznyakov 1975; Lyne & Graham-Smith 
1998). More recently, Percival et al. (1993) used ground- 
based optical broadband photometry together with a NUV 
photometric point from HST/HSP to determine a slope of 
a v = 0.11 ± 0.13, while Nasuti et al. (1996) obtained a 
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spectrum in the limited wavelength range 4900 - 7000 A 
(see below) and determined a v — —0.10 ± 0.01. 

Applying R = 3.1 and E(B-V) = 0.52 from §3.2 to our 
UV spectrum gives a spectral index of a v — — 0.035±0.040, 
where the error is simply the RMS around the fit. Obvi- 
ously, the main uncertainty in this procedure is the extinc- 
tion parameters. For E(B—V) in the range [0.48, 0.55], the 
uncertainty due to reddening becomes a„ — — 0.0351q'°3 4 . 
An accurate estimate of the spectral index obtained in the 
UV, for the given extinction range, can be obtained by 
a v = -0.035 + 3.12[E(B-V) -0.52]. Note that the spec- 
tral index for E(B — V) = 0.55 was erroneously reported in 
G98. The analysis in this paper supersedes this previous 
report. 

Including the optical spectrum from the NOT gives a 
wider wavelength range for the fit. The dereddened spec- 
trum shown in Figure 4 was obtained with R = 3.1 and 
E(B — V) = 0.52. The best power law fit to the complete 
spectrum (excluding Lya) gives a v = 0.11. 

Using the full wavelength range we have also tried to 
constrain both E{B — V) and R. This can be done us- 
ing the extinction curves from Fitzpatrick (1999), which 
are a one-parameter family in R. By assuming an intrin- 
sic power law we thus allowed both R and E(B — V) to 
vary, and chose the values that minimized o~ av , the stan- 
dard deviation of the power law fit. This procedure gives 
R = 3.0 and E(B — V) = 0.51, which is consistent with the 
values used above. In Figure 6 we show the dereddened 
spectrum of the Crab pulsar for several different values of 
R and E(B — V). In this figure we have also included IR 
data from Eikenberry et al. (1997). These plots clearly 
show the ambiguity of the reddening correction. 

We can express the spectral index of the power law 
fit for the values E[B - V) = 0.52 and R = 3.1 as 
a v = 0.11 ± 0.04±°;22- The first error represents the RMS 
around the power law fit, and the last errors include all 
power law fits in the extinction intervals R = [2.9,3.3], 
E(B-V) = [0.48,0.55]. A linear fit, a v = 0.11 - 0.38(i?- 
3.1) + 3.88[E(B - V) - 0.52] reproduces the obtained a v 
in this interval to better than 0.02 units. 

We can only echo Oke (1969) in his conclusion that the 
uncertainties in reddening corrections are large enough 
to allow both for negative and positive slopes. The in- 
ferred energy distribution for the electrons is given by 
7 = 0.8 ± 0.5. It is worth noting that the only other young 
pulsar for which a spectrum has been obt ained in the opt i- 
cal, PSR B0540-69, had a v = -1.6 ± 0.4 (|Hill et al. 19971) . 
This is clearly steeper than the Crab pulsar spectrum. 

Finally, the time tag mode of our FUV observations also 
allows us to extract phase-resolved spectra, and we have 
looked for spectral differences during the pulse phase. We 
found no significant differences above the 5% level, neither 
in the spectrum of the primary versus the secondary peak, 
nor in the leading versus trailing part of the primary peak. 

3.4. The 5900 A feature 

In the first report on the optical spectrum of the Crab 
pulsar, Minkowski (1942) reported a featureless continuum 
with no absorption or emission lines. The observations of 
Oke (1969) indicated a similar conclusion, although the 
spectral resolution was not very good. Since then, few 
attempts have been made to obtain a better optical spec- 



trum of the Crab pulsar, despite that technical develop- 
ment certainly admits improvements. The only modern 
optical spectrum of the Crab pulsar was obtained with 
the New Technology Telescope (NTT) by Nasuti et al. 
(1996). This phase-averaged spectrum covered a rather 
small wavelength region (4900-7000 A) and was flat with 
a v = -0.10 ± 0.01 if dereddened with E{B - V) = 0.51. 
The spectrum was also reported to show a large dip at 
~ 5900 A, the width being ~ 100 A. According to the au- 
thors this feature probably originates close to the pulsar 
itself, but no detailed physical mechanism was proposed. 

We have searched for this feature in our optical spectra, 
but found no sign for a dip. This is true for all the different 
gratings covering this region during all of the nights of data 
obtained in December 1998 (see Fig. 4). Neither do we see 
the feature in unpublished data taken by us at NOT one 
year earlier (December 1997). Nasuti et al. (1996) noted 
that the feature became enhanced after flux calibration. 
This suggests that the dip may be an artefact of the re- 
duction procedure. While their investigation used a single 
spectrophotometric standard star with flux sampled only 
every ioo A, we have used two standard stars with 2 A 
sampling in the relevant wavelength region. Although the 
feature could be time dependent, we propose to regard it 
as an artefact until confirmed by other observations. 

3.5. Absorption lines 

As in the NUV spectrum of the pulsar (G98), the FUV 
spectrum contains several absorption lines. The identified 
lines are shown in Figures 7 and 8, and their equivalent 
widths (EWs) and corresponding column densities, N, are 
listed in Table 3. For completeness, we have also included 
the lines in the NUV spectrum identified and measured by 
G98. 

The strongest line is Lya, and its damping wings can 
be used to estimate A(H I). To do this we assume that 
the optical depth in the damping wings is defined by 
ct(A)A(H I) (Shull & Van Steenberg, 1985). Here the ab- 
sorption cross section cr(A) = (4.26 x 10~ 20 cm 2 )/(A-A ) 2 , 
where the wavelengths are in A and A = 1215.67 A. 
We have not considered the instrumental profile, which 
is much narrower (~ 1 A) than the line width (see Fig. 7). 
From this analysis we obtain N(H I)= (3.0 ± 0.5) x 10 21 
cm" 2 . The spectrum corrected for Lya absorption is 
shown in Figure 7. Our estimated uncertainty in A(H I) 
reflects the uncertainty in the continuum fit of the cor- 
rected spectrum. The column density of free electrons 
toward the Crab is measured from pulsa r dispersion to 



be A c = (0.1755 ± 0.0007) x 10 21 cm" 2 (Cornelia et al 



1969). The amount of atomic hydrogen is thus jV(H)= 



(3.2 ± 0.5) x 10 21 cm" 2 . This value agrees well with both 
the result of Schattcnburg & Canizares (1986) who ob- 
tained N(B) = (3.45 ±0.42) x 10 21 cm" 2 from an estimate 
of the X-ray absorption toward the Crab, and the value 
we estimate from E(B — V) using the relation in de Boer, 
Jura, & Shull (1987), which gives N(R I)~ 3.0 x 10 21 cm" 2 
for E(B -V) = 0.52. 

The other absorption lines we identify around zero veloc- 
ity in the FUV spectrum are: C I AA 1277,1329,1561,1658, 
C II A1335, O I A1303, Al II A1671, Si II AA 1260,1527 
and Si IV A1394. The Si IV line is close to noise level, 
which explains why the Si IV A1403 component is not 
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seen. The C I AA 1277,1329 lines are also marginal de- 
tections, but their strengths correspond to those expected 
when compared with the strengths of C I AA 1561,1658. 
The C IV AA 1548,1551 doublet is absent at zero velocity, 
but shows a blueshifted absorption with a maximum shift 
of ~ 2500 km s _1 . We return to this line in §3.7. 

All lines except Lya, and the C IV doublet which is 
not interstellar or from a slow wind (§3.7), are unresolved 
by the moderate spectral resolution in the STIS spectra 
(w 1.2 A in FUV and w 3.2 A in NUV, corresponding to 
250 km s _1 and 400 km s _1 for the central wavelengths 
of the two gratings, respectively). The measured EWs of 
the absorption lines are rather large, which means that 
we cannot assume that the lines are resolved and optically 
thin (i.e., we cannot use "weak-line" theory [e.g., Morton 
1991]) to derive abundances of the various species, as this 
would severely underestimate the abundances. 

In the absence of a proper model for the distribution of 
intervening matter, we assume that the absorption is dom- 
inated by a single cloud component. We set the intrinsic 
"Doppler" width in this component to 1 km s _1 , and as- 
sume that this is the same for all lines. Furthermore, we 
consider only one spectral component of the lines listed in 
Table 3. All these assumptions cause us to systematically 
overestimate the column densities and abundances so that 
our estimates will be upper limits to these. With this in 
mind, we calculate Voigt line profiles and EWs as func- 
tions of column density, using the atomic data in Morton 
(1991). From the measured EWs we then obtain column 
densities and abundances for the interstellar gas toward 
the Crab. The abundances in Table 3 are presented in the 
standard logarithmic form where the value for hydrogen is 
set to 12.0. 

In Table 4, where we have simply coadded the abun- 
dances of the different ionization stages of each element, 
we give the abundances of C, O, Mg, Al, Si and Fe from our 
analysis. We also list the solar values for these elements 
according to Anders & Grevesse (1989). There appears to 
be no extreme depletion of any element, contrary to what 
was stated in G98, where "weak-line" theory was used. 
This is in agreement with the X-ray observations by Schat- 
tenburg & Canizares (1986) which were consistent with a 
solar abundance of oxygen. Our oxygen abundance from 
the single O I A1303 line is, however, rather uncertain due 
to geocoronal airglow corrections, and a possible blend- 
ing with Si II A1304. As pointed out above, our method is 
likely to systematically overestimate abundances. Broader 
lines, or more cloud components, would lower the abun- 
dances derived in accordance with depletion seen in the 
normal interstellar medium (e.g., 0.65 dex for carbon at 
2 kpc, Welty et al. 1999). Spectra with higher resolution 
are required to refine this analysis. 

3.6. Emission lines 

The FUV-spectrum contains also information about the 
Crab nebula. The long slit covers ~ 25" x 0'.'5 of the neb- 
ula, along position angle (PA) 40°. 5. The orientation of 
the slit is shown in Figure 1. To extract the pulsar spec- 
trum we used only 13 of the 1024 pixels of the MAMA 
detector. To obtain a spectrum of the nebular emission 
outside this extraction window we again used CALSTIS to 
produce rectified, wavelength- and flux-calibrated images. 
As the count rates were very low, we summed the emis- 



sion from two 8'.'2 long regions positioned 0'.'66 above and 
below the pulsar position. This safely excludes any contri- 
bution from the point spread function (PSF) wings of the 
pulsar. All emission from all 6 observations was combined 
and averaged to increase the signal-to-noise. We present 
this nebular spectrum in Figure 9, which shows only the 
wavelength region 1400 — 1700 A. This is to exclude the 
geocoronal lines of Lya and O I AA 1303,1356. The spec- 
trum is the average of the regions above and below the 
pulsar. For the spectrum extracted above (North-East of) 
the pulsar, the continuum level was flattened artificially. 
This is because we observed a rising continuum in the red 
part of this spectrum, which is probably due to contam- 
ination of the nearby star seen in Figure 1. Moreover, 
the count rates in the continuum are only ~ 7.5 x 10 -5 
counts pixel -1 s" 1 . At such low count rates, trends in the 
dark currents might influence the continuum. No dark im- 
ages are subtracted from the FUV-MAMA observations, 
because the dark images are known to be variable and to 
have low count statistics. Therefore, we will not discuss 
further the slope of the continuum emission from the neb- 
ula. 

Only two emission lines intrinsic to the Crab nebula can 
be seen, C IV A1550 and He II A1640. The measured inten- 
sities are - (2.1 ± 0.8) x 10~ 15 and - (1.4 ± 0.8) x 10~ 15 
ergs s" 1 cm~ 2 , respectively. This is an average for the 
two 8'.'2 x 0'.'5 regions, just above and below the pulsar. 
Our detection of these two lines is in accordance with the 
two previous spectroscopic observations of the Crab neb- 
ula in the FUV. Davidson et al. (1982) used the WE to 
detect these two lines, as well as C III] A1909. Blair et al. 
(1992) recalibrated the WE data and complemented these 
with observations from the Hopkins Ultraviolet Telescope 
(HUT). 

The widths of the lines are approximately 13 A and 
10 A, for C IV A1550 and He II A1640, respectively. The 
width is due to the velocity distribution of the material, 
the doublet nature of C IV A1550 as well as to the spec- 
tral resolution for extended sources for the spectrograph. 
The bulk of the emission is redshifted by ~ 1000 km s _1 , 
although a narrower zero- velocity component is seen from 
the region above the pulsar. 

The flux ratio of C IV A1550 to He II A1640 is thus 
~ 1.5, which is in accordance with the findings of David- 
son et al. (1982), who sampled the fluxes over a larger 
field, 20" x 10". They used this to argue that the Crab 
nebula has no overabundance of carbon. However, the 
subsequent investigation by Blair et al. (1992) showed 
that variations in the He II/C IV line ratio exist in the 
nebula. They cautioned on conclusions regarding the car- 
bon abundance, as differences in physical parameters (e.g., 
ionization, density, temperature, clumping) are difficult to 
disentangle from abundance variations. To determine the 
carbon abundance is important because it holds the poten- 
tial to provide information about the ZAMS mass of the 
progenitor. Comparisons to detailed photoionization mod- 
els for the large apertures used by WE and HUT are ham- 
pered by the fact that ionization con ditions arc known t o 
change over very small spatial scales ( (Bankrit et al. 1998 ) . 
The observations we present have finer spatial resolution 
and are not biased toward bright filaments. HST obser- 
vations of the same regions in the optical, to establish the 
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ionization conditions, could make more quantitative esti- 
mates of the carbon abundance possible. 

3.7. The outer shell 

As was mentioned in §3.5, and shown in Figure 8, 
C IV A1550 is the only line in the UV spectrum which 
shows clear evidence of blueshifted absorption. No evi- 
dence for absorption can be seen at zero velocity. The 
greatest absorption appears to arise in material moving 
at ~ 1200 km s^ 1 , and can thus be due to "normal" 
Crab nebula material. However, the absorption seems 
to continue to higher velocities, and can be traced out 
to ~ 2500 km s _1 (see Fig. 10). To estimate the sig- 
nificance of the detection we have calculated an average 
and RMS from two 50 A regions redward and blueward 
of the C IV line. Although the line is too noisy for indi- 
vidual pixels to be significant, it consists of many pixels 
that are consistently below the average. The eleven pixels 
between 1650 — 2780 km s" 1 give a detected absorption 
with a significance of 5.8<r. For the five points in the region 
2330 — 2780 km s _1 the significance is 3. Oct. 

Together with Clark et al. (1983), who also reported 
velocities in excess of 2000 km s _1 , this is the highest ve- 
locity ever measured in the Crab, and can be interpreted 
as evidence for the existence of the long sought fast outer 
shell (cf. §1). The column density of C IV in Table 4, 
N(C IV) = (3.0 ±1.1) x 10 14 cm" 2 , assumes the line to be 
resolved and optically thin. This should be the case for ve- 
locity broadening caused by the tentative fast shell, though 
we caution that the material moving at <~ 1200 km s -1 
may not be spectrally resolved. The total column density 
of C IV is therefore likely to be higher, while that of the 
proposed fast shell is lower than the value given in Table 
4. We will now investigate whether the detection is con- 
sistent with a fast outer shell model and what information 
can be provided about the supernova ejecta. 

3.7.1. Constraints from C IV X1550 

We adopt a model similar to that of LFC86, i.e., outside 
the observed filaments we attach a massive, freely coasting, 
spherically symmetric shell. The inner radius, i?; n , is set 
to 5.0 x 10 18 cm, which agrees with the "mean" inner ra- 
dius of the presumed shell used by Sankrit & Hester (1997, 
see their Fig. 7). For free expansion this corresponds to a 
velocity of w 1680 km s _1 . Furthermore, we assume that 
the mass of the shell, M s h, is 4 Mq, and that the outer 
radius of the shell, R out , is 1-9 x 10 19 cm. The maximum 
velocity, V out , is then « 6370 km s _1 , and the total kinetic 
energy of the shell, E s h, is 1.0 x 10 51 ergs, if the density is 
constant in the shell. For a density which decreases with 
radius, E s ^ is lower, if M s h is held constant. This is shown 
in Table 5 for density slopes up to 77 = 9, where r] is defined 
as p(R) = p(R in ) (R/R in )-'i. 

In our model we have used the relative abundances 
X(K): X(Rc): X(C) = 1.0: 0.1: 3.5 x 10~ 4 , where the 
X(C)/X(H) ratio corresponds to the solar value. With 
these abundances we have calculated the absorption in 
C IV AA 1548,1551 as a function of n and the relative frac- 
tion of carbon in C IV, X(C IV). The parameter X(C IV) 
is unity when carbon is all in C IV. 

In Figure 10 we show results for two models where 
X(C IV) has been kept constant throughout the shell. The 
dotted line shows the C IV line in the case of 77 = and 



X(C IV) = 1, while the dashed line is for 77 = 3 and 
X(C IV) = 0.14. Both models are described in Table 5. 
The photoionization models of LFC86 show a case simi- 
lar to our 77 = case (their model 1). There, carbon is 
ionized beyond C IV, and only a very small fraction, con- 
centrated to the outer edge of the shell, remains in C IV. 
This gives maximum absorption in the 77 = case around 
-6000 km s" 1 (see Fig. 1 of LFC86 for such a case). This 
is clearly not what is seen in the FUV data. On the con- 
trary, the observed absorption peaks at lower velocities, 
meaning that the C IV number density must be the high- 
est closer to i?; n . The photoionization models of LFC86 
make the 77 = case rather unlikely from this point of 
view. Instead we turn our attention to steeper density 
profiles. 

Such a situation is highlighted by the 77 = 3 case in Fig- 
ure 10. Our assumption of constant X(C IV) in the shell 
is probably more realistic for 77 = 3 than for 77 = 0. This 
is because the ionization parameter, £ = n 7 /n , where n 7 
and n e are number densities of ionizing photons and elec- 
trons, respectively, has the radial dependence £ cx R v ~ 2 , if 
absorption can be neglected. With absorption included, 
77 must be > 2 to obtain a near-constant £,(R). Our 
choice of X(C IV) = 0.14 has only been made to fit the 
data. Any combination of X(C IV) and M s h, so that 
X(C IV)M sh w 0.56 Mq, would fit the data equally well. 
In our model, for 77 = 3, this gives a minimum mass (for 
X(C IV) = 1) of the shell of M sh - 0.6 Mq, and a mini- 
mum kinetic energy of E s h ~8x 10 ergs. 

While a constant X(C IV) in the shell gives a good fit 
to the line profile for 77 = 3, cases with a steeper den- 
sity dependence must have an ionization structure where 
X(C IV) increases outward through the shell. From 
the radial dependence on the ionization parameter this 
could be possible, but may require a fine-tuning between 
the R n ~ 2 and exp(— r) parts of £. Time dependent ef- 
fects are also likely to be important, and the likelihood 
of an increasing X(C IV) with radius can only be ex- 
plored by numerical models. We postpone this to a future 
study. To give the same absorption at 2500 km s _1 as the 
77 = 3 model in Figure 10, models with other 77 must have 
X(C IV) ~ 1.49 1 '- 3 X(C IV) in . The parameter X(C IV) in 
is the value of X(C IV) at i?i n , and is the value given in 
Table 5. For example, for 77 = 9, X(C IV) must be ~ 11 
times higher at the radius of 2500 km s _1 compared to 
V(CIV) in . 

It thus appears as if the observed line profile of the 
C IV A1550 absorption is consistent with the outer shell 
being fast and massive. That is, the C IV line cannot ex- 
clude that the shell could carry an energy of 10 51 ergs, be- 
cause the absorption in the gas with the highest velocities 
(~ 6000 km s _1 ) in the outer shell model would simply 
disappear in the noise of our rather poor signal-to-noise 
spectrum. 

To limit ourselves to the observed velocities, we adopt 
a maximum velocity of 3000 km s _1 to obtain lower limits 
on the mass and energy of the shell. Table 5 shows that 
M s h could be very low (i.e., M s h = M trunc X(C IV), with 
M tmnc defined in Table 5) if 77 is large and if we only use the 
absorption at i?j n as a criterion. However, if we also require 
that the absorption at 2500 km s _1 should be the same 
for any 77 as in the 77 = 3 model shown in Figure 10, the 
minimum mass for models with 77 > 3 is given by M s h ;> 
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lAW' 3 X(C IV) in M trunc Mq, which for 77 = 3 (4, 5, 7, 
9) becomes M sh £ 0.27 (0.31, 0.37, 0.55, 0.90) M . The 
kinetic energy corresponding to this M s ^ (i.e., for V^, ut = 
3000 km s" 1 ) is E sh £ 1.5 (1.7, 1.8, 2.4, 3.4) x 10 49 ergs. 
The limits on M s h and -E s h for shallower density profiles 
are fixed by the product X(C IV)M trunc (cf. Table 5), and 
are M sh £ 0.48 (0.36) M & and E sh £ 3.1 (2.2) x 10 49 ergs, 
respectively, for r\ = 1 (2). All these limits scale inversely 
with the overall abundance of carbon. 

To summarize the constraints from C IV A1550, we first 
emphasize that the line shows that an outer shell with 
maximum velocity of V^ ut ~ 2500 km s _1 appears to be 
present. We have a ;> ha detection of velocities between 
1650 — 2780 km s _1 . The kinetic energy of the shell de- 
pends on its mass, density structure and extent. For a rel- 
atively shallow density profile (p oc i?~ 3 , or shallower) the 
kinetic energy of the shell could be as high as 10 51 ergs, if 
the shell is spherically symmetric and extends to velocities 
higher than those we can detect with the obtained signal- 
to-noise. For a maximum velocity of ~ 6000 km s~ 4 , the 
mass of the shell required to obtain this kinetic energy 
would be 4 — 8 Mq. A completely flat density profile, 
however, seems unlikely from the results of LFC86, both 
in the fast case (LFC86, model 1), and for a model similar 
to our V out w 3000 km s" 1 case (LFC86, model 4). 

For a density profile steeper than p oc i?~ 3 the kinetic 
energy is most likely < 10 51 ergs since the shell mass 
should be lower than ~ 8 Mq to be compatible with pro- 
genitor models. The lowest mass and energy we estimate 
for a shell with V ou t = 3000 km s _1 is for 77 ~ 3, and 
are ~ 0.3 Mq and ~ 1.5 x 10 ergs, respectively. These 
values are approximate as they depend on spherical sym- 
metry, the value of R{ n , and a model fit to the line profile 
of C IV A1548 at -2500 km s" 1 where the line profile 
is rather uncertain. To distinguish between models with 
different density slopes (see Table 5), photoionization cal- 
culations are needed. 

3.7.2. Other constraints on the shell 

There are, unfortunately, no other lines detected in the 
UV data that can constrain our current analysis further. 
Two potentially useful doublets are Si IV AA 1394,1403 
and N V AA 1239,1243, but we cannot identify absorp- 
tion at high velocities in any of these two doublets. The 
absence of the two doublets is, however, not surprising. 
Silicon should be more highly ionized than carbon, and 
LFC86 found that C IV A1550 should produce significantly 
stronger absorption than N V A1240. The spectral region 
around N V A 1240 is also rather noisy and the line sits 
in the damping wing of Lya. Although the absence of the 
lines cannot constrain models in this simple analysis, it can 
be used in conjunction with photoionization calculations 
to test different models. 

The tentative outer shell has been previously searched 
for also in the optical. Searches in [O III] have been neg- 
ative for the region outside the observed [O III] skin (cf. 
Fesen et al. 1997, and references therein). This is not 
surprising from the point of view of the models of LFC86, 
where oxygen is more highly ionized than O III. 

A highly ionized massive shell is bound to give rise 
to Ha emission. The deepest search for such emission 
was done by Fesen et al. (1997) who found a surface 



brightness limit in Ha of 1.5 x 10 7 ergs cm 2 s 1 sr 1 . 
With the dereddening suggested by Fesen et al. (1997), 
i.e., Ahq = 2.536E(B - V), and E(B - V) = 0.52 (cf. 
above), the dereddened surface brightness limit becomes 
5.1 x 10~ 7 ergs cm~ 2 s _1 sr -1 . We have calculated the 
surface brightness in our models to see how it compares 
with the observed limit. We use a temperature of the ten- 
tative shell of 2 x 10 4 K, which is nearly three times higher 
than that used by Fesen et al. (1997) and Murdin (1994), 
but in accordance with the models of LFC86. 

The maximum (dereddened) surface brightness in Ha, 
Eh q , occurs at the impact parameter, p — R ln , i.e., just at 
the edge of the observed nebula. The value of £h~ q f° r this 
impact parameter is given as a function of 77 in Table 5 for 
the model with M s h = 4 M Q . It is seen that the modeled 
surface brightness exceeds the observed limit for 77 ;> 4. 
However, Shq decreases rapidly with p for large 77. Table 
5 shows that at p = 1.1 x R m (corresponding to ~ 17" 
outside the observed nebula for a distance of 2 kpc, and 
roughly where the search by Fesen et al. (1997) was con- 
ducted), the surface brightness exceeds the observed limit 
only for r\ J> 5. For a shell mass as low as 0.3 — 0.9 Mq, 
which we found to be likely lower limits to the shell mass 
for the density slopes investigated, the shell would easily 
have escaped detection in Ha. 

A method to derive parameters for the outer shell was 
devised by Sankrit & Hester (1997). They estimated the 
density needed to form a radiative shock at the interface 
between the nebula and the presumed outer shell, as such 
a shock is needed in their model to explain the observed 
[O III] skin. They estimate that a minimum density of 
yo/mn ~ 12 cm -3 is needed, at least in the presumed equa- 
torial plane of the nebula. If this is true also in the direc- 
tion toward the pulsar, the models of Sankrit & Hester 
yield N(C IV) - 10 14 cm" 2 for the radiative tail of the 
shock. This would escape detection in our data since the 
intrinsic line width in their model should be small (much 
less than our spectral resolution). We can therefore not 
distinguish between a radiative or adiabatic shock (or no 
shock at all) in the direction to the pulsar. This also means 
that it is unlikely that any of the absorption we detect oc- 
curs in a region similar to the radiative region in the model 
of Sankrit & Hester. 

Sankrit & Hester (1997) estimate that a mean density 
at the inner edge of the shell, averaged over all polar 
angles, should be ~ 8 cm~ 3 . This would correspond to 
"■H~(^in) ~ 5.7 cm~ 3 for the He/H ratio used in Table 5. 
To see if we can make a consistent model including the 
C IV line, the Ha surface brightness limit by Fesen et al. 
(1997) and the model by Sankrit & Hester (1997), we have 
assumed an upper limit to the shell mass of 8 Mq (inside 
V = 3000 km s _1 ), and used the information in Table 5. 
We then find that n H (i?in) ^ 5.7 (M trunc /8 Mq) cm" 3 is 
required to get a high enough density at R la to agree with 
the model of Sankrit & Hester (1997). According to the 
values in Table 5 this is fulfilled for 77 ;> 3. The mass and 
kinetic energy for a shell with such a high 7Jh (R[ n ) and with 
V out = 3000 km s" 1 would be 6.8 (5.0, 3.8, 2.4, 1.7) M Q and 
3.9 (2.7,1.8,1.0,0.7) x 10 50 ergs, respectively, for 77 = 3 
(4, 5, 7, 9). A caveat for this model is that Sh q at 
p = 1.1 x R ln then becomes 2.0 (1.4,1.0,0.60,0.35) x 
10~ 6 ergs cm~ 2 s _1 sr -1 for 77 = 3 (4, 5, 7, 9), which 
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is close to, or higher than, the observed limit. In this 
scenario, a limit on the Ha surface brightness improved 
by a factor of a few, close to the observed nebula, should 
be able to distinguish between models with different den- 
sity slopes even if the He/H ratio is higher than we have 
assumed. Photoionization models to accurately calculate 
the temperature and to check the radial dependence on 
X(C IV), are also needed. 

4. CONCLUSIONS 

Using STIS onboard the HST we have observed the 
Crab nebula and its pulsar in the far-UV (1140 - 1720 A). 
We have obtained the pulse profile of the pulsar, which 
is very similar to our previous near-UV profile, although 
the primary peak appears to be marginally narrower than 
in the near-UV data (5%, 2a). Combining the far- and 
near-UV data, and assuming an intrinsic power law for 
the pulsar continuum, we have derived an extinction of 
E(B — V) = 0.52 mag toward the Crab. No evidence for 
a non-standard extinction curve was found. We have also 
added optical spectra taken with the NOT to obtain a 
spectrum of the pulsar from 1140 A to 9250 A. We have 
shown that the pulsar spectrum can be well fitted over 
the full UV/optical range by a power law with spectral 
index ct v = 0.11. The exact value of the spectral index 
is, however, sensitive to the amount and characteristics of 
the interstellar reddening, and we have investigated this 
dependence for a likely range of E(B — V) and R. In 
the optical, we find no evidence for the dip in the pulsar 
spectrum around 5900 A reported by Nasuti et al. (1996). 

The interstellar absorption lines detected in the UV have 
been analyzed, and are consistent with normal interstel- 
lar abundances. The column density of neutral hydrogen 
is (3.0 ± 0.5) x 10 21 cm" 2 , which corresponds well to the 
value derived for E(B — V). From the Crab nebula itself 
we detect the emission lines C IV A1550 and He II A1640. 
The ratio of the fluxes of these lines is similar to what 
has been derived previously, although obtained with much 
improved spatial resolution. 

C IV A1550 is also seen in absorption toward the pul- 
sar. The line is broad and blueshiftcd with a maximum 
velocity of ~ 2500 km s _1 , and there is no absorption at 
zero velocity. These are the highest velocities measured in 
the Crab and shows that there exists material outside the 
visible nebula. This can be interpreted as evidence for the 
fast shell that has been predicted to surround the Crab 
nebula (Chevalier 1977). We have used a simple, spher- 
ically symmetric model in which the density in the shell 
falls off with radius as R~ v from 5 x 10 18 cm (correspond- 
ing to w 1680 km s" 1 ) to derive the mass and energy of 
such a shell. The conclusions from our model depend on 
how we tie our model into other observations and mod- 
els. From the C IV line alone, we find that the minimum 
mass and kinetic energy of the fast gas are ~ 0.3 M and 
~ 1.5 x 10 ergs, respectively. This occurs for a density 
slope r] ~ 3. A model with a flat (ry = 0) density profile 
appears unlikely as the required ionization structure dis- 
agrees with the modeling of Lundqvist et al. (1986). The 
maximum mass of the shell is set by progenitor models, 
and is unlikely to be much larger than 8 M . With a high 
shell mass, and the velocity extending to velocities much 
higher than we can detect, the shell could carry an energy 
of 10 51 ergs. The signal-to-noise of C IV A1550 is too low 



at high velocities to reject or confirm such a conclusion. 

Adding constraints from the model of Sankrit & Hester 
(1997) to those from the C IV line narrows down the pa- 
rameter space for the shell. In particular, a density slope 
of rj ;> 3 is required to agree with the interpretation of the 
observed [O III]-skin being a radiative shock. For rj < 9, 
the shell mass is then ^ 1.7 M and the kinetic energy 
J> 7 x 10 49 ergs. Although the limit on the Ha surface 
brightness from the search of Fesen et al. (1997) tend 
to favor models with steep density profiles, a model with 
r/ = 3 might still be possible, if the He/H ratio is higher 
than solar also in the fast shell, and the asymmetry of the 
outer shell different from that in the model of Sankrit & 
Hester (1997). 

We thank Rob Fesen for help during preparations of 
the HST observations, and for discussions and comments 
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Table 1 

LOG OF STIS FUV OBSERVATIONS 



Observation 


Start time (MJD) 


Exposure time 




(51200.0+) 


(seconds) 


O4ZP01010 


.54924 


2100 


O4ZP01020 


.60619 


2460 


O4ZP01030 


.67384 


2460 


O4ZP01040 


.74103 


2460 


O4ZP02010 


.81802 


2100 


O4ZP02020 


.87488 


2460 



Table 2 

LOG OF NOT OBSERVATIONS 



Grism 


Date 


Wavelength range 


Resolution 




(Doc 1998) 


(A) 


(A) 


3 


24+25 


3300-6400 


6.6 


5 


24+25 


5200-9250 


9.0 


6 


27 


3200-5500 


4.7 


7 


25 


3820-6100 


5.0 


8 


25 


5830-8340 


3.6 



Table 3 
ABSORPTION LINES 



Species 


A a 


f a 

Jose 


Equivalent width 


Log (Column density) 

(cm- 2 ) 


Log (Abundance) 




(A) 




(A) 




H I 


1215.67 


0.416 


37 ±5 


21.48 ± 0.08 


12.0 


C I 


1277.46 


0.0966 


0.18 + 0.12 


17.49 ± 0.71 


7.96 ± 0.76 


C I 


1329.34 


0.0580 


0.11 + 0.11 


< 18.11 


< 8.71 


C I 


1561.05 


0.0804 


0.35 + 0.17 


18^20 ± 0.46 


8.70 + 0.48 


C I 


1657.59 


0.140 


0.54 + 0.17 


17.89 + 0.29 


8.40 ± 0.30 


C II 


1335.31 


0.128 


0.46 + 0.24 


18.15 ± 0.50 


8.65 ± 0.53 


C IV 


1549.05 


0.286 


1.84 + 0.65 


14.45 ± 0.16 b 


N/A b 


O I c 


1303.49 


0.0488 


0.54 + 0.21 


18.51 ± 0.36 


9.02 ± 0.38 


Mg I 


2852.96 


1.830 


0.93 ± 0.33 


16.11 ± 0.34 


6.62 ± 0.35 


Mg II 


2796.35 


0.612 


0.79 ± 0.25 


16.79 + 0.28 


7.30 ± 0.30 


Mg II 


2803.53 


0.305 


0.55 + 0.27 


16.70 + 0.47 


7.20 ± 0.49 


Al II 


1670.79 


1.833 


0.56 ± 0.34 


16.04 ± 0.57 


6.54 + 0.60 


Si II 


1260.42 


1.007 


0.60 + 0.18 


16.82 ± 0.28 


7.33 + 0.29 


Si II 


1526.71 


0.230 


0.32 + 0.14 


16.57 + 0.39 


7.08 ± 0.41 


Si IV 


1393.75 


0.514 


0.33 + 0.20 


16.63 ± 0.63 


7.12 + 0.66 


Fe II 


2586.65 


0.0646 


0.27 + 0.20 


17.15 + 0.86 


7.58 + 0.96 


Fe II 


2600.17 


0.224 


0.40 + 0.21 


16.71 ± 0.53 


7.21 + 0.56 



a Morton (1991) 

b Not interstellar. The column density is calculated from "weak-line" theory (Morton 1991). 

c Could be affected by subtraction of gcocoronal emission, and is probably also blended with Si II A1304.37 
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Table 4 
ABUNDANCES 



Element Solar value a This paper b 



C 8.55 8.9 ±0.4 

O 8.87 9.0 ±0.4 

Mg 7.58 7.4 ±0.3 

Al 6.47 6.5 ±0.6 

Si 7.55 7.6 ±0.4 

Fe 7.51 7.2 ± 0.6 



a Anders & Grevessc (1989) 

b Our method is likely to systematically 
overestimate the abundances. See §3.5. 



Table 5 

PROPERTIES OF A 4 M n FAST SHELL. a 



„b 


n H (Rin) b 




Eh. 


C 


S Ha d 


X(C IV) C 


Mtrunc 




(0.714 X p/m H , cm -3 ) 


(10 51 orgs) 


(ergs cm 2 


s- 1 sr- 1 ) 


(ergs cm 2 s 1 sr 1 ) 




( M ) 





0.12 


0.99 (0.03) s 


9.3(4.1)« 


x 10~ 9 


9.2(3.9) g X 10~ 9 


l h 


0.43 


1 


0.32 


0.86 (0.05) 


2.4(1.8) 


x 10~ 8 


2.1(1.6) X 10~ 8 


0.62 


0.78 


2 


0.77 


0.72 (0.08) 


8.1(7.6) 


x 10~ 8 


6.1(5.6) x 10~ 8 


0.28 


1.29 


3 


1.62 


0.56 (0.11) 


2.7(2.7) 


x 10~ 7 


1.7(1.6) x 10~ 7 


0.14 


1.92 


4 


2.94 


0.42 (0.14) 


7.4(7.4) 


x 10~ 7 


3.8(3.8) x 10~ 7 


0.081 


2.57 


5 


4.66 


0.32 (0.15) 


1.6(1.6) 


x 10~ 6 


6.9(6.9) x 10~ 7 


0.053 


3.11 


7 


8.71 


0.21 (0.16) 


4.7(4.7) 


x 10~ 6 


1.4(1.4) x 10~ 6 


0.030 


3.71 


9 


13.0 


0.17 (0.15) 


9.1(9.1) x 10~ 6 


1.8(1.8) x 10~ 6 


0.021 


3.92 



a Minimum and maximum radii arc 5 X 10 18 cm and 1.9 X 10 19 cm, respectively, corresponding to the velocities 1680 km s 1 
and f» 6370 km s — 1 . 

b n H (R) = n u (R in ) (R/R in )-v 

c Ho surface brightness at impact parameter p = i?; n , i.e., just at the edge of the observed nebula. (T = 2 X 10 4 K). 

d Ha surface brightness at impact parameter p = l.lij; n , i.e., ~ 17" outside the observed nebula. (T = 2 X 10 4 K). 

e Fraction of carbon in C IV to obtain the observed optical depth in C IV A1550 at ~ 1680 km s _1 . (X(C)/X(H) = 3.5 X 10" 4 ). 

'How much of the 4 Mq shell that resides inside 3000 km s _1 . Density and X(C IV) at i?j n are assumed to be the same as for 
the 4 Mq shell with maximum velocity 6370 km s~ - 1 . 

g Valucs in parentheses are for the case with the shell only reaching 3000 km s _1 , i.e., for the masses in the last column in the 
table. 

h For r] = the maximum optical depth in C IV A1550 is too small even with X(C IV) = 1. (See Fig. 10). 
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Fig. 1. — An R-band image of the Crab nebula obtained at the NOT in December, 1998. Shown are the extents and position angles of the 
slits used in the optical and in the FUV. The STIS-FUV slit is 25" long. The NUV observation (G98) used a 2"x2"aperture centered on the 
pulsar. 
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Fig. 2.— Pulse profile of the Crab pulsar in the FUV (full line). Also shown is the NUV pulse profile from G98 (dashed line). The FUV 
and NUV data were processed in the same way. The blow-up of the primary peak shows the FUV to be slightly narrower. The statistical 
errors range from 0.003 in the valleys to 0.014 at the peak for both the FUV and NUV profiles, as illustrated by the (±1<t) error bars in the 
insert. Note that the count rate of the pulsar signal is close to zero in the interpulse region. 
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Fig. 3. — Spectrum of the Crab pulsar in the UV. This is a combination of our FUV (1140 — 1730 A) observations obtained in January 
1999, and the NUV (1600 — 3200 A) data from G98. The spectra overlap nicely. The Lya absorption dip shows some residuals from the 
subtraction of strong geocoronal emission. 
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Fig. 4. — Spectrum of the Crab pulsar in the UV and in the optical. The optical data are from NOT. The lower spectrum shows the 
flux-calibrated spectrum without dereddening. The optical spectrum connects rather well to the NUV data. The uppermost spectrum has 
been dereddened with E(B — V) = 0.52 and R = 3.1. The full line shows the best power law fit, which has the spectral index a v = 0.11. The 
insert shows a blowup of the optical region where Nasuti et al. (1996) reported a broad absorption dip in the spectrum. 
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Fig. 5. — The UV extinction curve in the direction toward the Crab. This curve was obtained by assuming the pulsar intrinsic spectrum 
to follow a power law. For comparison, the mean galactic extinction curve from Fitzpatrick (1999) is also shown. 
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Fig. 6. — The Crab pulsar spectrum dcreddened with different values of E(B — V) and R. The upper panel shows the dereddened spectrum 
for E(B — V) = 0.52 and three different values of R: 2.9 (upper), 3.1 and 3.3 (lower). The 1R data points are from Eikenberry et al. (1997) 
and are not included in the power law fits. The fitted power laws have a v = 0.19, 0.11 and 0.041, respectively. The lower panel shows 
dereddening for R = 3.1 and E(B — V) =0.49 (upper), 0.52 and 0.55 (lower). The fitted power laws have spectral indices a„ = —-0.005,0.11 
and 0.23, respectively. The individual spectra in both panels have been shifted by —0.3, 0.0, and +0.3 dcx in the vertical direction for clarity. 
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Fig. 7. — FUV spectrum around Lya. The dashed line shows the original data, and the solid line shows the spectrum after compensation 
for absorption in the damping wings of hya. The column density of neutral hydrogen toward the Crab determined in this way is N(R I) = 
(3.0 ± 0.5) X 10 21 cm~ 2 . The straight dashed line shows a continuum fit to the spectrum after correction for Lya absorption. 
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Fig. 8. — Absorption lines seen in the spectrum toward the Crab pulsar. Only C IV A1550 shows absorption which cannot be interstellar 
(see Fig. 10). The difference in velocity between the two components of the C IV doublet is shown by the two vertical lines marked by "AA". 
The equivalent widths of all lines are given in Table 3. Note that the zero velocity of the C II and C IV multiplets in the figure corresponds 
to the most blueward component of the multiplets, while Table 3 lists weighted wavelengths. 
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Fig. 9. — FUV emission from the Crab nebula in the region around the pulsar. The spectrum was obtained by averaging the emission 
from two 8'.' 2 long region above and below the pulsar position along the 0"5 slit. The continuum has been flattened out and the spectrum 
smoothed with a 3 pixel boxcar average. 
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Fig. 10. — Pulsar spectrum around C IV A1550. No absorption at zero velocity of C IV A1550 is seen. Instead the line is blueshifted, 
reaching a maximum velocity ~ 2500 km s _1 . The spectrum also shows the interstellar lines Si II A1527 and C I A1561, neither of which 
blends with the C IV line. Si II A1533, which sometimes accompanies Si II A1527 in interstellar spectra, appears to be absent, and cannot 
explain the observed absorption at velocities < 2700 km s — 1 . Overlaid on the observed spectrum are two of the models for the fast shell in 

Table 5: r\ = (dotted) and r\ = 3 (dashed), where r\ is defined from p(R) = p{R ln ) {R/Rin)~ rl > an d ^in is the inner radius of the shell. We 
have assumed that R; n = 5 X 10 18 cm, which corresponds to a minimum shell velocity of *=s 1680 km s _1 . See Table 5 for details of these and 
other models of the shell. 



